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Abstract 

We discuss the Type-X (lepton-specific) two Higgs doublet model as a solution 
of the anomaly of the muon g — 2. We consider various experimental constraints on 
the parameter space such as direct searches for extra Higgs bosons at the LEP H 
and the LHC Run-I, electroweak precision observables, the decay of Bg —)• pA, 
and the leptonic decay of the tau lepton. We find that the measurement of the 
tau decay provides the most important constraint, which excludes the parameter 
region that can explain the muon g — 2 anomaly at the lo" level. We then discuss 
the phenomenology of extra Higgs bosons and the standard model-like Higgs boson 
(h) to probe the scenario favored by the g — 2 data at the collider experiments. We 
find that the 4r, 3r and At + W/Z signatures are expected as the main signal of 
the extra Higgs bosons at the LHC. In addition, we clarify that the value of the 
hTT coupling is predicted to be the standard model value times about — 1.6 to — 1 . 0 , 
and the branching fraction of the /i —)■ 77 mode deviates from the standard model 
prediction by —30% to —15%. Furthermore, we find that the exotic decay mode, 
h decaying into the Z boson and a light CP-odd scalar boson, is allowed, and its 
branching fraction can be a few percent. These deviations in the property of h will 
be tested by the precision measurements at future collider experiments. 



1 Introduction 


The anomalous magnetic moment of the muon = {g — 2)/2, so-called muon g — 2, 

is a very precisely measured observable. The latest measurement of by the E821 

collaboration [1] gives 

a;fP = 11 659 208.0 (5.4)(3.3) x 10"^°. (1) 

As it has been well known that there is a discrepancy between the experimental value 
and the prediction of the standard model (SM). According to the calculation evaluated 
in Refs. [2,3] 

= (28.7 ± 8.0) x 10“^°, (Davier et. ah) 

= (26.1 ± 8.0) X 10“^°, (Hagiwara et. al.) 

the discrepancy is more than the 3a level, which can be considered as an indirect evidence 
of the existence of a new physics model. This discrepancy will be further probed at 
Fermilab [4] and J-PARC [5] in the near future. Since the size of the deviation is the 
same order as the electroweak contribution = 15.4 x 10“^° [6], we expect that new 
physics exists at the electroweak scale if the strength of new interactions is as large as 
that of the weak interaction. In such a new physics scenario, new particles are expected 
to be light enough to be directly discovered at the LHC. Therefore, it is quite interesting 
to consider models beyond the SM as a solution of the muon g — 2 anomaly. 

Among various models which can explain the anomaly (for a review, e.g., see Ref. [7]), 
two Higgs doublet models (2HDMs) give simple solutions. In 2HDMs, there are extra 
Higgs bosons {H, A, and H^) in addition to the SM-like Higgs boson [h), and they can 
give new contributions to a^. Usually, a softly-broken discrete Z 2 symmetry is imposed [8] 
to avoid flavor changing neutral current (FCNC) processes at the tree level. Under the Z 2 
symmetry, four independent types of Yukawa interactions are allowed depending on the 
assignment of the Z 2 charge to the SM fermion [9,10], which are called as Type-I, Type- 
H, Type-X (or lepton specihc) and Type-Y (or flipped) [11]. In all the types of Yukawa 
interactions, the lepton couplings to the extra Higgs bosons can be sizable enough to 
explain a^. In the Type-I and Type-Y 2HDMs, however, the top Yukawa coupling also 
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becomes large together with the enhancement of the lepton couplings. This is disfavored 
from the view point of perturbativity. Thus, the Type-II and Type-X 2HDMs are suitable 
to solve the muon g — 2 anomaly. 

The muon g — 2 has been calculated in a number of papers within 2HDMs [12-22]. 
In the early 2000s, this was calculated at the one-loop level in the Type-II 2HDM in 
Ref. [12]. After that, it was pointed out in Refs. [13,14] that the two-loop Barr-Zee type 
diagrams [23, 24] give a signihcant contribution to if a mass of A is (9(10-100) GeV 
and if there is large Abb or At^t~ couplings. In Ref. [19], the implication on collider 
signatures was discussed in the Type-X 2HDM, namely, the h —)■ AA —)■ 4r process can 
be important in the favored parameter region by After the discovery of the Higgs 
boson at the LHC [25,26], the muon g — 2 has been reanalyzed by taking into account 
the Higgs boson search data in addition to the previous experimental constraints [20-22]. 
Furthermore, the recent observation of Br(R 5 —)■ gAg,~) at the LHC [27] gives a new 
constraint on the parameter space of 2HDMs [21]. 

The difference between the Type-II and Type-X 2HDMs is the quark couplings to 
the extra Higgs bosons. In the Type-II 2HDM, both the lepton and down-type quark 
couplings are enhanced simultaneously, and thus the model is severely constrained by 
flavor physics and direct searches of the extra Higgs bosons. On the other hand, in the 
Type-X 2HDM, the quark couplings to the extra Higgs bosons are suppressed when the 
lepton couplings are enhanced. Thus, the constraints are weaker than those in the Type-II 
2HDM. In fact, it was clarihed in Refs. [20,21] that only the Type-X 2HDM can solve the 
muon g — 2 anomaly with satisfying the current experimental data^. 

Another important constraint comes from the lepton flavor physics. In the Type-X 
2HDM, the constraint from the leptonic r decay [11,30-32] gives a severe constraint on 
the parameter space favored to explain the g — 2 anomaly because of the lepton coupling 
enhancements. However, this important constraint has not been included in the previous 
analyses. Therefore, in this paper, we calculate the leptonic r decay and the Ztt vertex 
at the one-loop level in the Type-X 2HDM in order to compare the precise experimental 
measurements. We then investigate the favored parameter region by under these 

^In addition to the muon g — 2 anomaly, there are several other motivations for this model. For 
example, see Refs. [28,29]. 
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Table 1: Assignment of the Z 2 parity and factors in Eq. (22) in each type of the Yukawa 
interactions. 

constraints in addition to those already known. Furthermore, we evaluate the running 
of the scalar quartic couplings by renormalization group equations (RGEs), and require 
that the couplings do not become too large up to a certain energy scale, for example 10 
TeV. We hud that extra Higgs boson loop contributions can reduce the discrepancy in 
to be 2(7 level, but not less than la level. We then study the collider phenomenology in 
the favored parameter region. 

This paper is organized as follows. In Sec. 2, we define the Lagrangian of the 2HDM, 
and derive the Higgs boson couplings with the gauge bosons and the fermions. In Sec. 3, 
we discuss constraints from direct searches for the extra Higgs bosons at LEP H and the 
LHC Run-I, electroweak precision observables, the decay of —)■ , the leptonic r 

decay, and the triviality bound. In Sec. 4, we show the favored parameter regions by the 
muon g — 2 anomaly. In Sec. 5, we discuss the collider phenomenology of the extra Higgs 
bosons at the LHC, the deviations in the SM-like Higgs boson h couplings, and the decay 
branching fractions of h. We also discuss the exotic decay mode: h ^ ZA. Conclusion 
is given in Sec. 6 . In Appendix, we present the expressions for the decay rates of extra 
Higgs bosons and those for the parton level cross sections for the production of extra 
Higgs bosons at the LHC. 

2 The 2HDMs 

In this section, we dehne the Lagrangian of the 2HDM, in which the Higgs sector is 
composed of two SU{2)l doublet scalar fields Hi and H 2 . To avoid the tree level FCNC, 
we impose a Z 2 symmetry in the Higgs sector which can be softly-broken in general. Under 
the Z 2 parity, four types of Yukawa interactions are dehned depending on the assignment 
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of Z 2 charge as listed in Table 1. 

The most general Higgs potential with the softly-broken parity is given as 




(m^, 2 ^lH 2 + h.c.) 


Ai 


A 


+ -^\H2\^ + \s\H,\^\H2\^ + \,\HIH2\^ + 


+ fc,c. 


( 2 ) 


Throughout the paper, we consider the CP-conserving case of the Higgs sector for sim¬ 
plicity, so that the imaginary parts of 771^2 and A 5 are assumed to be zero. The Higgs 
helds are parametrized as 


H^ = 




K 

+ hi 


itti 


(* = 1 , 2 ), 


(3) 


where Vi and V 2 are the VEVs of the Higgs doublets which are related to the Fermi 
constant Gp by + v 2 = 1/{\/2Gf)- The ratio of the two VEVs is parametrized 

by tan [5 = V 2 /V 1 . 

The mass eigenstates of the scalar bosons are expressed by introducing the mixing 
angles a and (3 as 


hi 

h2 


ai 

0-2 


hf 

ht 


cos a — sin a \ / H \ 
sin a cos a ) \ h^ ) ' 

cos/d — siu/d A / A 
siu/d cos/d y y H y ’ 

cos /d — sin /3 A / G^ A 
siny cos/3 


(4) 

(5) 

( 6 ) 


where G^ and G^ are the Nambu-Goldstone bosons which are absorbed by the Z and W 
bosons as the longitudinal component, respectively. 

The squared masses for the physical Higgs bosons are given by 

m^± = 

m\ = — v^As, 

mjj = cos^(q; — + sin^(Q; — /d)M 22 -|- sin2(Q; — /3)M^^2, 

ml = sin^(Q; — + cos^(q; — /d)M |2 — sin2(Q; — /d)Mf 2 , (7) 
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where = 771 ^ 2 /(sin cos /3) describes the breaking scale of the Z 2 symmetry, and 
are given by 

v‘^ 

^11 = n^(Ai cos"^ /d + A 2 sin"^ f3) + y A 345 sin^ 2/3, ( 8 ) 

M 22 = + n^(Ai + A 2 — 2 A 345 ) sin^ /d cos^ /d, (9) 

M ^2 = y (A 2 sin^/d — Ai cos^/d + A 345 cos2/d) sin2/d, (10) 

where A 345 = A 3 + A 4 + A 5 . The mixing angle a is also expressed in terms of Mf^ as 

2 M^ 

tan 2 («-/^) = 

All the quartic coupling constants in the Higgs potential can be rewritten in terms of the 
physical parameters as 

Aiu^ = —tan^ [3 + (m^ tan^ /3 + + (m^ + ml tan^ I3)c^p-a 

+ 2{m]j - ml) tan (3sf}-aCf}-a, 

\2v‘^ = cot^ 13 + (m|- cot^ [3 + ml)sl_^ + {m^j + ml cot^ I3)cl_^ 

- 2{mjj - ml) cot f3sis_c,C0-a, 

Xsv^ = 2m]j± - + {ml - - (tan/3 - cot 13) 

+ m\ — 2 m^±, 

X^v^ = — m\, ( 12 ) 

where s_s_q, = sin(/d — a) and cp-a = cos(/d — a). 

A size of some combinations of A’s in the Higgs potential is constrained by taking into 
account perturbative unitarity [33-36] and vacuum stability [37,38]. Through Eq. (12), 
such a constraint can be translated into a bound on the physical parameters; e.g., the 
masses of the scalar bosons. First, the condition for vacuum stability; he., the requirement 
for bounded from below in any direction of the Higgs potential with large scalar fields, is 
given by [37,38] 

Ai >0, A 2 > 0, -s/ A 1 A 2 + A 3 + MIN(0, A 4 + A 5 , A 4 — A 5 ) > 0. (13) 

Second, the perturbative unitarity bound is obtained by requiring that all the eigenvalues 
of the s-wave amplitude matrix for the elastic scatterings of two body boson states 
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are satisfied as 


K±\ < 2- (14) 

All the independent eigenvalues were derived in Refs. [34-36] as 

~ 3(Ai + A 2 ) ± \/9(Ai — A2)^ + 4(2A3 + A4)^ , (15) 

3z7r L J 

Ogi = (-^1 + -^ 2 ) ± (Ai — A2)^ + 4A| , (16) 

®3,± = (-^1 + -^ 2 ) ± Y^^(Ai"-Ov2)^”+4^ , (17) 

4± = ^(A3 + 2A4±3A5), (18) 

<± = J^(V±A4 ), (19) 

<± = ( 20 ) 

The Yukawa interaction terms are given by 

>^^Yukawa = “ VuH'^Qlu'r - ydHlQidR - yeHlLie^ + h.c., (21) 

where Hu = H^Hu- In Eq. (21), Hu, Hd and Hi are either Hi or H 2 depending on the 
type of Yukawa interaction. In the mass eigenstates of the scalar bosons, the interaction 
terms are expressed as 


-'Yukawa 




f=u,d,i 


V2K,H*u (d - 

V V J 


H+uPni + h.c. , (22) 


where T| = +1/2 (—1/2) for f = u {d,i), and Vffi is the Cabibbo-Kobayashi-Maskawa 
matrix element. The and factors are dehned by 

'^y—a T ^fC/d—a, Cy—a Cf-^y—a- (^3) 

The factors in Eqs. (22) and (23) are given in Table 1. 
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Figure 1: One-loop (left) and two-loop Barr-Zee (right) diagrams which give corrections 
to the muon g — 2. 


From the kinetic terms of the scalar helds, the ratios of the coupling constant among 
the CP-even scalars and gauge bosons are extracted as 


ghvv 

ghVV,SM 


Qhvv 

Cf^—a- 

ghVV,SM 


{V = W, Z) 


(24) 


As it is seen in Eqs. (22), (23) and (24), in the limit of sin(/d — a) —?• 1, both hff and 
hVV couplings become the same as those in the SM, so that we can call this limit as the 
SM-like limit. 


3 Constraints on the Type-X 2HDM 

In the 2HDMs, the one-loop diagrams and the two-loop Barr-Zee type diagrams shown 
in Fig. 1 give dominant contributions to the muon g — 2. It has been known that the 
Barr-Zee type diagrams give a sizable positive contribution to in the case of a large 
coupling and a small as pointed it out in Refs. [13,14]. In the Type-X 2HDMs, 
a large At^^~ can be realized by taking tan /d 3> 1 since = — tan f3 as shown in Table 1. 
Typically, when tan/? > 40 and rriA = 0(10-100) GeV, the muon g — 2 anomaly can be 
explained in the Type-X 2HDM [20]. In this section, we focus on the Type-X 2HDM 
with the large tan /? and small mA scenario to explain the g — 2 anomaly, and we discuss 
important experimental constraints in this situation. 
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3.1 Direct searches for the extra Higgs bosons 

There has been no signal of the extra Higgs bosons at any collider experiments. This 
gives lower limits on the masses of the extra Higgs bosons depending on the magnitude 
of couplings with SM particles. We hrst summarize the current bounds from the LEP H 
experiment, and we also review those from the LHC Run-I. 

3.1.1 LEP II 

There are constraints on the masses of the extra Higgs bosons from the direct production 
at the LEP H experiment with the maximal collision energy to be about 200 GeV. From 
the pair production process e’''e“ —)■ —)■ H^H~ the lower bound was obtained 

by mH± > 93.5 GeV at 95 % C.L. [39] under the assumption of Br(iL+ —)■ T^Vr) = 1 
which is realized by tan /3 >2, sin(/3 — a) ~ 1 in the Type-X 2HDM. 

From the pair production of the neutral Higgs bosons e’'“e“ —)■ Z* —)■ AH, the lower 
bound for the sum of uia and uih is given to be about 190-195 GeV for niA > 40 GeV [40] 
under the assumption of Br{H —)■ r+r“) = Br(H —)■ = 1 which is realized by 

tan /3 >3 and sin(/d — a) ~ 1 in the Type-X 2HDM. 

The searches for A and H from the bremsstrahlung process e’''e“ —)■ ttA/H have also 
been performed for the range of rriA/H < 50 GeV. This process gives an upper bound on 
tan (3 for a fixed value of niA/H- For example, tan [3 > 128.1 (120.1) and tan (3 > 44.8 (40.0) 
are respectively excluded at 95% C.L. for {jnu) to be 30 GeV and 15 GeV [41] with 
the case of Br(H —)■ rr)=Br(iL —)■ rr)=l. 

We note that the branching fractions for the extra Higgs bosons into a fermion pair 
can be reduced when there is a non-zero mass splitting among them. For example, —)■ 
HWb) and H — )■ open in the case of mH± > "nr a and mn > ttia, respectively. There 

also happen the inverse processes like A —)■ iL=*=H/^b) A —)■ HZ^*^ as long as they are 
kinematically allowed. In such decay modes associated with a gauge boson, the bounds 
on masses on the extra Higgs bosons can be weaker than those given in the above. 
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3.1.2 LHC Run-I 


At the LHC, extra Higgs boson searches have been performed in various channels. In 
the most of channels, an enhancement of the Yukawa couplings of the extra Higgs bosons 
becomes important to obtain a bound on their masses or coupling constants. However, 
in the Type-X 2HDM, the couplings of the neutral extra Higgs bosons to the quarks are 
suppressed by cot (3. Thus, the processes such as gg ^ A/H ^ rr and gg —)■ bbA/bbH —)■ 
bbTT do not set a limit on the masses in a large tan f3 case. 

Similar to the neutral Higgs boson productions, the cross section of the production 
such as gb —)■ H^t is also suppressed by cot^ /3 in the Type-X 2HDM. If mH±+mb < rrit, the 
top decay t —)■ H^b can be used to constrain mH±- From the process pp ^ tt ^ bbH^W^ 
with —)■ r^z/, the upper limit on BR(f —)■ H^b)xBIi{H^ —)■ r^z/) has been driven to 

be between 0.23% and 1.3% at 95% C.L. for mH± in the range of 80 GeV to 160 GeV [42]. 
This gives the bounds, for example, tan/3 < 6 and 15 for mH± = 100 and 150 GeV at 


95% C.L. in the Type-X 2HDM using 0.23% of the product of the branching fractions. 

Apart from the production processes via Yukawa couplings, one must take care of the 
h —)■ AA decay in the case of tua < In the Type-X 2HDM, this typically gives the 

four r hnal state, because the A ^ tt decay can be the main decay mode as explained 
in Sec. 3.1.1. In Ref. [43], the upper bound on Br(h —)■ AA —)■ 4r) is given to be about 
0.2 for rriA > 30 GeV and 0.2-0.5 for 15 < tua < 30 GeV. In the 2HDMs, the branching 
fraction is determined by the dimensionless hAA coupling XhAA dehned as the coefficient 
of the hAA vertex in the Lagrangian; he., C = vX^aa hAA + ■ ■ ■ which is given by 


1 


XhAA = TT^ [(2M^ - 2m\ - ml)sp_a -f (M^ - m^)(cot/3 - tan/3)c;3_„ 


2u2 

The partial decay width of h —)■ AA is then expressed by 

.,2 


h^AA — 


\2 

■^hAA^ 

Snnih 


^ Fsm X 


m 


h 


XhAA \ ' 

0.015 / 


4m^ 


m 


2 ’ 


(25) 


(26) 


where Tsm = 4.41 MeV is the total decay width of the SM Higgs boson for = 
125 GeV [44]. Therefore, to satisfy Br(h —)■ AA)< 0.2, XhAA^^A x 10“^ is required. 
We can simply take XhAA = 0 by setting an appropriate value oi (3 — a from Eq. (25) as 


tan(/3 


M 2 


ml 


a) = 


2M2 


2m\ 


mi 


-(tan (3 — cot (3). 


(27) 
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In the case of tan/3 3> 1, m\lm‘\j± 1 and m\/ra\j± 1, we obtain 


sin(/3 — a) ~ 1 — 
2 

cos(/3 — a) ~ 


tan^ /3 


mf 


tan /3 


1 + 


1 + 

m 

2m]j± 


2m\ 


H± 


m 


H± 


m: 


m 


H± 


( 28 ) 

(29) 


From the above expressions, we hnd that the SM-like behavior of h is realized by taking 
tan /3 3> 1, because of sin(/3 — a) ~ 1. 


3.2 Electroweak precision observables 

The extra Higgs bosons can modify the electroweak precision observables from the SM 
prediction via the loop effects. Such an effect can be used as an indirect search for the extra 
Higgs bosons and also used to constrain parameter space in the 2HDM. In this subsection, 
we discuss the constraints from the oblique parameters and the Z boson decay. 


3.2.1 Oblique parameters 


The electroweak oblique S, T and U parameters are introduced by Peskin and Takeuchi [45] 
which parametrize new physics effects on the gauge boson two point functions. These pa¬ 
rameters are calculated in 2HDMs in Refs. [46-52]. In the SM-like limit sin(/3 — a) —)■ 1, 
these parameters are given to be the same formulae as those given in the inert doublet 
model [53]. For the case of uia mz <^mH± — rnn, the contribution from the additional 
scalar bosons is given by 


AT ~ 


- ~ 0 . 022 , 

7271 

mH{mH±-mH) _ f mn \ ( mH± - mH \ 

327r2aemt^2 - V200 Gevy V 10 GeV ) ' 


(30) 

(31) 


We also hnd that AH is the same order as AS* in our setup for large AT regime. If we 
take tuh = and sin(/3 — a) = 1, the Higgs potential respects the custodial SU{2)v 
symmetry [48,50], which makes AT = 0. The S and T parameters driven by the Ghtter 
group [54] are 


A^ = 0.05 ±0.11, AT = 0.09 ±0.13, 


(32) 
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with the reference values of mh = 125 GeV and rrit = 173 GeV. The prediction of AS 
parameter is inside the la error of the measured value, and the T parameter constrains 
on the mass splitting \mH — \ = O{10) GeV. Hence we take uih = rnH± to avoid the 

constraint on the oblique parameters in the following analysis. 

3.2.2 Z boson decay 

The property of the Z boson such as the mass, the total width, and the decay branching 
ratios were precisely measured at the LEP experiment. If new physics particles modify 
such a precisely measured quantity, their masses and/or couplings are severely constrained. 

In our scenario, the Zt~^t~ vertex can be signihcantly deviated from the SM prediction 
by loop effects of the extra Higgs bosons, because they strongly interact with charged 
leptons in the large tan/? case. In order to discuss how the modihed vertex affects the 
observables, we dehne the effective Zff vertex as 

C = gzlYivf - a/75)/Z^, (33) 

where gz = ol cosO^ and 6^ being the weak mixing angle. Although there are several 
dehnitions for sin^ 9w, we here use the on-shell dehnition [55] of it which is determined by 
using W and Z boson masses, he., sin^ 6 'vk = 1 — rn\rlm\. The effective vector coupling 
Vf and axial vector coupling df can be separately written by the contributions from the 
tree level and from the one-loop level as 

■0/= u /-|-df = af + af°^, (34) 

where the tree level contributions are expressed as 

Vf = - sin^ OwQf, af = ^Tf, (35) 

with Qf being the electric charge of /. The loop contributions and are composed 
of the counter term and the one particle irreducible (IPI) Zff vertex correction diagram: 

Uy °^ = -f ^loop _ - 1 - ^ipi^ 
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After imposing the on-shell renormalization conditions, the counter term contribution is 
expressed by [56] 


Svf = - Vf 


[n//y(p^) + ^ff;siP^)]p^=m} 


Self — — (If 


d 


n//,u("^/) + [n//,y(P ) + ^ff,siP )]p 2 =m 2 




[m 


fh 


(37) 

(38) 


where diagram contributions to the fermion two point functions dehned 


as 


uff\p^) = mffW) - ( 39 ) 

In the SM-like limit sin(/3 — a) = 1, the deviation in and purely comes from the 
extra Higgs boson loop diagrams. In this case for / = r, we obtain 


Au]°°P = - u]°°P(SM) 


16n^ 




VrFi{mH) + VrFi{niA) - 2arFi{mH±) + {v^ + a^)F2(m//±,m^i) 


(40) 


A4°°p = a[°°P - a[°°P(SM) 


IGtt^ V V 


—6 


a^Fi(mH) + a^iPi(mA) - 2 a^Fi(m//±) 

+ {vr + ar)F2{niH±,niH±) - AarF2{niH,niA) 


where the loop functions are given as 

,2 K 1 7"! 


^ 5 1 

FAm) = ln^ + --m + - 


m 


2 ■ 4 -■■ ■ 2 


dx 


'0 


F 2 (mi,m 2 ) = ^(ln^ + l) - 
4 V mf 


r-l 


dx 


d ~ 1) + ^^(2 + y — yx) 

Jq m|a;(l — y) — m?{l — x) 

y'^ — m|(l — 2|/)] 


(41) 

(42) 

(43) 


n / JO Jo ‘2 ml{l - y) + xy[mj - ml{l - y)] 

In the above expressions, we neglect the mass of the tau lepton in the loop functions. We 
note that the ^ 2 ( 7711 , m 2 ) function is invariant under the interchange of mi <H- m 2 , so that 
Au]°°P and Aa[°°P does not change the value under uih t-)- niA- 

Let us apply the modihed vertex to the leptonic partial decay width of the Z 

boson Z —)■ ii: 


T{Z^ii)=gl^{v‘i + d^,). 


(44) 
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We define the ratio of the partial width of Z —)■ rr to that of Z —)■ ee as 

= r(Z TT)/r(Z ee). (45) 

The deviation in the ratio from the SM predictions are then given by 

ARr/e = Rr/e “ ^ [t'£ReAnJ°°P + aiReAa°°^] . (46) 

' ' ^ r(Z-)■ ££)sm Stt l ^ j 

The SM prediction is given by 

T{Z ii)sM = 83.995 ± 0.010 MeV. (47) 

The measured values of the leptonic decay width and Rr/e are given by [60] 

r(Z ^ ££)exp = 83.984 ±0.086 MeV, = 1.0019 ± 0.0032. (48) 

We hnd that tan /3 > 50 (70) is excluded for rriA = 10 (50) GeV when uih = mH± = 
300 GeV. The bound becomes weaker for the larger mH± ■ We will combine the constraint 
from the Z ^ ii decay with the muon g — 2 result in Sec. 4. 

3.3 Flavor experiments 

Effects of the extra Higgs bosons can appear in various observables measured at flavor 
experiments. Therefore, similar to the electroweak precision measurements, flavor mea¬ 
surements can be used to constrain the parameter space in the 2HDM. In this subsection, 
we discuss Bg —)■ and the leptonic decay of r. 

3 . 3.1 Bg —)■ jjA jj, 

It was pointed out in Ref. [57] that the branching fraction of the Bg —)■ process in 
the Type-II 2HDM is enhanced by the factor tan"^ (5 which comes from the box type and 
the penguin type diagrams with extra Higgs boson mediation. The origin of the tan^/d 
dependence is that both the charged lepton and down-type quark couplings to the extra 
Higgs bosons are proportional to tan in the Type-H 2HDM. On the other hand, in the 
Type-X 2HDM, the lepton couplings are enhanced by tan 13, while the quark couplings 
are suppressed by cot (3. Thus, the tan^ (3 dependence does not appear in the branching 
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fraction of Bg —)■ and the additional leading contribution is almost independent on 

tan f5 for large tan 13. Although the deviation from the SM becomes mild as compared to 
the case of Type-II, we check Bt{Bs — )■ because the light CP-odd Higgs boson A 

could give a sizable contribution, which is required to explain muon g — 2 anomaly. 

For the calculation of Bg —)■ gAgT in the 2HDM, we use the formulae given in Ref. [58]. 
The observed branching fraction at the LHC is Br(Rs —)■ g~) = (2.9 ± 0.7) x 10“® [27]. 

We show the constraint from Bg —)■ g^g~ on the parameter space of 2HDM in Sec. 4. 


3.3.2 Leptonic r decay at the tree level 


In the SM, the leptonic r decay is caused by the W boson exchange diagram at tree level. 
In the 2HDM, the mediated diagram also contributes to the leptonic r decay. The 
effect of contribution on the partial decay width of r was calculated in Refs. [11,31], 
and that on the Michel parameters, which is dehned just below, in Ref. [32]. 

The differential decay rate of r —)■ gv^Vr is given in terms of the Michel parameters 
{p,ri,S and and G^r dehned in Eqs. (52) and (53) as [59] 


where u = (m^ -|- m‘j^)/2mr, x = and xq = m^/u with being the muon energy. 

Vr is the polarization of the tan, and 6 is the angle between the polarization and the 
momentum direction of the muon. The functions F{x) and A{x) are dehned as 


2 

F{x) = x{l — x) -I—p(4a;^ — 3x — x^) -|- 77 x 0(1 — x), 
9 

A{x) = 1 — X + ^5(4x — 4+ \/l — xl). 

O * 


By using x = mgrriT-taiA (3 we hnd^ 


P = 


3 

4’ 


/ ;^2x1/2 

G,r = GFl^l + jj , 

2x 3 

4 + ^2’ “4’ ^ “ 4 + ■ 


^We find these expressions are inconsistent with Ref. [32] 


(50) 

(51) 


(52) 

(53) 
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Z = m, / Z = ^^ m, / 

Figure 2: The ratio of the decay rate (upper left and upper right) and the Michel param¬ 
eters r] (lower left) and ^ (lower right) as a function of 2 ; in the leptonic tan decay. In 
each panel, outside regions from the dotted lines are excluded at the 95 % C.L. 
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We see that p and 5 are equal to the SM values at the tree level. The observed Michel 
parameters of the r decay are p = 0.013 ± 0.020 and ^ = 0.985 ± 0.030 [60]. The ratio of 
the decay rate in the 2HDM to that in the SM prediction is given as [11,31] 


_ r(r —)■ 1 _ 2 

Gf ) ~ r(r /iz/z/)sM ^TUrf {ml/ml) 4’ 


(54) 


where the phase functions f{x) and g{x) are given hy f{x) = 1 — 8a; — 12a;^ log x + 8x^ — x'^ 
and 5 f(a;) = l + 9a; —9a;^ —a;^ + 6a;(l + a;) logo;. To hnd a constraint on Eq. (54), we can use 
the constraint on the flavor universality. In the similar manner to Eq. (54), we introduce 
Gefi and Ger- Since nie, nif^ m-r, the corresponding terms to the rightest term in Eq. (54) 
for GefM and Ger are 1, and thus Ge^ = Ger = Gj? in 2HDM. There are constraints on the 
lepton universality given by HEAG group [61]^ 

G G 

^ = 1.0029 ±0.0015, -^ = 1.0018 ±0.0014, (55) 

^er 


and their correlation coefficient is 0.48. We use this bound and Eq. (54) to make constraint 
on 2HDMh 

In Fig. 2, we show the dependence of the ratio of the decay rate given in Eq. (54) 
(upper two panels) and the Michel parameters p (lower left) and ^ (lower right). First, 
from the upper panels we can see that the allowed ranges of are found to be 2 : < 0.003 
and 0.50 ^ z < 0.57. Second, from the lower left panel, > 0.05 is excluded by the 
measurement of p. The constraints from the ^ parameter is weaker than that from p. 
Therefore, by combining the first and the second statements, the allowed region of is 
restricted to be 2 : < 0.003 By using ~ 1.88 x 10“^ x (tan/d/30)^ x (300 GeV/m^±)^, we 
hnd that tan /3 > 38 is excluded for mH± = 300 GeV. 


3.3.3 Lepton universality at the one-loop level 

As we discussed in Sec. 3.3.2, the typical size of the contribution to the ratio of the 

tan decay is G(10“^) at the tree level as it is seen in Fig. 2. However, the SM prediction 

^The ratio of the effective Fermi constant G^r/Ge^ and Ger/Ge^ are corresponds to g^/ge, 

gr/ge and gr/g^i in Ref. [61], respectively. 

^At the tree level discussion, we can replace Gg^ and G^t with Gp, then we find G^T-jGp = 1.0023 ± 
0.0012 or [G^tIGfY = 1-0046 ± 0.0025 by combining Eq. (55). However, this treatment does not work 
for the calculation including one-loop corrections which is given in Sec. 3.3.3. 


17 



r 



T 



Figure 3: The leading one-loop diagrams for the leptonic tan decay process. 

is given at almost the lower edge of the experimental bound, (see Eq. (55)), so that the 
negative contribution to G^r/Gp of order is constrained. Thus, we focus on the 
quantum corrections to the process via W exchange diagram. 

The dominant contribution arises from the diagrams with picking up two tau Yukawa 
couplings which are proportional to (m^/u^) tan^/3. We show the diagrams which give 
the dominant contributions to the process at the one-loop level in Fig. 3. Other diagrams, 
such as box diagrams, are smaller than these contributions and we ignore them in this 
analysis. The quantum correction is flavor dependent, and there is no flavor dependent 
counter terms in this model, so the correction is finite. We find the contributions from 
Fig. 3 modihes the W-t-p^ couplings 


9Wtu —t Qw TV (1 + ig) . 


(56) 


where 



(57) 


In niA ~ rriHO case 



(58) 


Finally, we find 



^1- xl{2{A{l + 5gf + z‘^) - 8 xo2:(1 + 5g) 


- 5xo(4(l SgY + z^) - l<oxlz{l + Sg)^ 
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mH± = 150Ge V m//± = 1 OOGe V 
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Figure 4: Constraints on the (m^, tan/9)-plane by the leptonic tan decay at the one-loop 
level for mH± =300 (upper left), 200 (upper right), 150 (lower left), and 100 GeV (lower 
right) in the case with mu = mH±- The green, yellow, and red regions are excluded at 
90%, 95%, and 99% C.L., respectively. 


where u and Xq are dehned in Sec. 3.3.2. Using Eqs. (54) and (59), we hnd that the large 
tan /3 region and light charged Higgs region are disfavored, see Fig. 4. This constraint is 
the most severe for the explanation of the muon g — 2 anomaly, as we will see in Sec. 4 

3.4 Triviality bound 

In order to avoid the constraints from the various observables, we need to take large mass 
differences between A and and A and H. As a result, the Higgs quartic couplings are 
as large as 0{1). Such a large coupling can be grown up in a certain energy scale, and it 
becomes too strong to rely on the perturbative calculation. We thus take into account the 
triviality bound in which we require that all the Higgs quartic couplings do not exceed a 
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Figure 5: The triviality bound with A = 10 TeV (left) and A = 100 TeV (right). The 
numbers shown in the panel are Amax = 47 r and a/^. 

certain value until a given energy scale. 

We calculate the /^-functions up to the two loop level for the RGE by using SARAH [62], 
and run the couplings to higher energies. We treat the coupling values at the tree level 
as the input parameters at p = rrit. We dehne Amax = max{|Ai|, IA 2 I, IA 3 I, IA 4 I, IA 5 I} and 
the cutoff scale of the model, A, as Amax(A) = dvr or V^vr. 

The result for A = 10 TeV with requiring Amax < 47 r is shown in the left panel in 
Fig. 5. We hnd that < 370 GeV for tua — 20 GeV is required for A 10 TeV. This 
constraint on rn^± is stronger than the one from the perturbative unitarity bound using 
Eqs. (14)-(20), he., mH± ^ 700 GeV. We check that this result is consistent with that 
given in [63]. If we require Amax < instead of Amax < dvr, then the bound becomes 
stronger. In such a case, we hnd mH± — 260 GeV is required. Here we take tan/3 = 30 
aX jj, = rrit as the input value, but the result is insensitive for tan (3. We also plot the case 
for A = 100 TeV in the right panel in Fig. 5. The bound is stronger than A = 10 TeV 
case. We hnd mH± ^ 310 GeV (240 GeV) for Amax = dvr (a/Ht). Our parameter choice 
here is the same as in Sec. 4. 
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4 Muon g — 2 


We show the numerical results for the muon g — 2 under all the constraints discussed 
in the previous section. We calculate the muon g — 2 hj using 2 HDMC 1.6.4 [64] which 
contains the one-loop diagrams [12] and the two-loop Barr-Zee diagrams [13] as shown 
in Fig. 1. In the Type-X 2HDM, the contributions from the one-loop diagrams and the 
Barr-Zee diagrams are comparable and have opposite sign. Thus, we have to take into 
account both contributions. The contributions from the beyond the SM part should be 
sizable to solve the muon g — 2 anomaly, and thus at least one of the new particles has to 
be light. The lower mass bound on is order of 100 GeV, so that it can not be arbitrary 
small. The effect from H gives a destructive contribution to the SM one, so that it makes 
the situation becomes even worse if H is light. On the other hand, the effect from A 
gives the constructive contribution, and it makes the discrepancy small. Therefore, A is 
required to be lighter than H and in order to solve the muon g — 2 anomaly. 

The Higgs sector in the 2HDM has eight parameters. Two of them are hxed to 
reproduce the SM parameters, he., Gp = 2~^/‘^v~‘^ = 1.166379 x 10“^ GeV“^ and rrih = 
125 GeV. To suppress Br(/i —)■ AA), we set Xhaa = 0. In addition, we take = mH± 
to avoid a large contribution to AT. Furthermore, we £x Ai = 0.1. This Ai value is 
realized by taking ~ rn^j in the large tan g case, and the fixing the value of Ai is not 
significant to the result in the following. Therefore, we have three remaining parameters 
which can be expressed as tan/3, ttia and mH±- We note that in this parametrization, 
sin(/3 — a) is given as the output parameter, which is determined via Eq. (25). As it is 
seen in Eq. (28), 1 — sin(/3 — a) is suppressed by 1/ tan^ /3 in the large tan (3 case, so that 
h behaves as the SM-like Higgs boson. 

In Fig. 6, we show the prediction of on the m^-tan/3 plane in each fixed value of 
m^±, e.g., m^± = 200 (upper-left), 250 (upper-right), 300 (lower-left) and 350 GeV (lower 
right). In the dark and light shaded regions, can be explained at the la and 2a levels, 
respectively. We find that the measurement of r —)■ gives the stringent constraint 

on the parameter space in the m^-tan /3 plane. This constraint is getting stronger in the 
cases of smaller mH± or and larger tan/3. For example, in the case of mH± = 200 
GeV, tan/3 > 35 is excluded at the 95% G.L. in the case of uia = 20 GeV. For ~ 10 
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muo = mn+ = 200 GeV my^o = mn+ = 250 GeV 




m^o = roy^* = 300 GeV mf _|0 = 111 ^+ = 350 GeV 




[GeV] [GeV] 


Figure 6: Results in the Type-X 2HDM in the case of Xhaa = 0 and Ai = 0.1. Dark 
and light blue shaded regions can explain the muon g — 2 anomaly [3] at the la and 2a 
levels, respectively. We take mH±{= mn) =200, 250, 300 and 350 GeV in the upper-left, 
upper-right, lower-left and lower-right panels, respectively. The left region from the red 
line is excluded by the measurement of Bg —)■ p/x. The above regions of green, black, 
purple line are excluded by the r decay, the direct search at the LEP and the Z —)■ rr 
decay, respectively. All of the exclusions are given at the 95% C.L. 
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GeV, Br(i?s —)■ gives additional excluded regions which are not excluded by the 

measurement of the tan decay. The tan (3 dependence in By{Bs —)■ is negligible as 

we expected in the discussion in Sec. 3.3.1. In most of the parameter region shown in this 
hgure, the deviation of Bt{Bs —)■ from the SM value is 0(1) %. The constrained 

from Z —)■ rr is weaker than that from the tan decay in all the parameter regions shown 
in this hgure. 

Consequently, the parameter region which can explain the g — 2 anomaly at the la 
level is excluded by the measurement of the tan decay at the 95% C.L., and at best we 
can explain the anomaly at the 2a level. The typical parameters to explain muon g — 2 
anomaly at the 2a level is 10 ^ niA ~ 30 GeV, 200 ^ ~ GeV and 30 ^ tan (3 < 50. 

In this parameter space, however, the region with mj{± ^ 270 GeV has tension with the 
signal strength for the h ^ tt mode measured at the LHC as we will see in the next 
section. 

5 Impact on the Higgs phenomenology at collider ex¬ 
periments 

In the previous section, we have seen that the relatively light extra Higgs bosons and large 
tan f3 are favored to explain the muon g — 2 anomaly. Such a light particle can be directly 
discovered at the LHC Run-H and the International Linear Collider (ILC). Furthermore, 
the precise measurement of the property of the SM-like Higgs boson h will give an indirect 
probe of this scenario. In this section, we hrst discuss the decay and production of the 
extra Higgs bosons at the LHC, and then we investigate how the property of the SM-like 
Higgs boson is modihed in the favored parameter set indicated by the muon g — 2 in the 
Type-X 2HDM. 

Throughout this section, we consider the case of 

200 GeV < mH± < 400 GeV, = M = mH±, 10 GeV < < 30 GeV, 

30<tan/d<50, tan{[3 — a) =- -^ (tan/? — cot/d). (60) 

2M^ — 2mi — mi 
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In addition, the following values of the SM parameters are used [60,65]: 

mz = 91.1876 GeV, mw = 80.385 GeV, Gp = 1.1663787 x lO"® GeV^ 

mt = 173.07 GeV, rrih = 3.0 GeV, nic = 0.677 GeV, nir = 1.77684 GeV. (61) 

TTic and mi, are evaluated at the mz scale in MS scheme, which are taken from Ref. [65]. 
Other parameters are taken from Ref. [60]. The mass observed Higgs boson is taken to 
be 125 GeV. 

5.1 Phenomenology of the extra Higgs bosons 

First, we discuss the branching fraction of the extra Higgs bosons in the parameter set 
given in Eq. (60). For the GP-odd Higgs boson A, only the A ^ tt and A ^ bb modes are 
allowed at the tree level. Since the decay rate of the former (latter) channel is enhanced 
(suppressed) by tan^/d (cot^/d) in the Type-X 2HDM, the branching fraction of H —)■ rr 
becomes almost 100% in our scenario. Similar enhancement happens in the decay rates of 
77 —)■ rr and —)■ r^z/. At the same time, the other modes H —)■ AZ and —)■ AW^ 

are also important because of the large mass difference between H /and A. Therefore, 
the following decay modes should be taken into account; 

r^z/, AW^, H^tt, H ^ AZ, A ^ rr. (62) 

The formulae for the decay rates are given in Appendix A. We here note that the 77 —)■ AA 
and H ^ hh decays also open whose decay rates are determined by the trilinear HAA 
and Hhh couplings, respectively. However, these couplings are proportional to cos(/d —a), 
and its magnitude is suppressed by cot /7 in the large tan [3 case. Thus, these decay modes 
are not important in our scenario. 

In Fig. 7, we show the contour plots for the branching fractions of 77 —)■ AZ (left 
panel) and —)■ AW^ on the mH±Aaji (3 plane. The H ^ rr and 77+ —)■ r+z/ modes 

fulhll almost the remaining branching ratios for 77 and 77+, respectively. In the large mH± 
and small tan /7 region, the branching fractions of 77 —)■ AZ and 77+ —)■ AIF+ are getting 
larger. For example, we obtain Br(77 —)■ AZ) > 70 (90)% in the case of tan/7 < 35 with 
mjj = 200 GeV (tan/7 < 33 with mjj^ = 300 GeV), and Br(77+ —)■ AIF+) > 70 (90)% in 
the case of tan/7 < 38 with m^j = 200 GeV (tan/7 < 35 with mjj± = 300 GeV). 


24 



50 


Br(H -AZ) 


Br(H^- AW"^) 




Figure 7: Contour plots for the branching fractions of if —)■ AZ (left panel) and —)■ 

AW^ (right panel) on the m//±-tan/3 plane. We take = mii± and = 20 GeV. 
Br(ii —)■ 74Z)+Br(ii —)• rr) ~ 1 and Br(ii’'' —)■ 2lhF’'')+Br(ii’'' —)■ ~ 1 are satished. 

Next, we discuss the production process of the extra Higgs bosons at the LHC. As we 
mentioned in Sec. 3.1.2, the production processes via the Yukawa interaction cannot be 
used in the large tan (3 case in the Type-X 2HDM. Therefore, the following electroweak 
processes give the dominant production mode: 

pp Z* I-i* ^ ii+if-, pp^Z* ^ HA, 

pp^W* ^ H^H, pp^W* ^ H^A. (63) 

The analytic expressions for the parton level cross section are given in Appendix B. We hnd 
that the cross sections are determined by the masses of the extra Higgs bosons and sin(/d — 
a), and they do not depend on the type of Yukawa interactions. By using CalcHEP [66] 
with CTEQ6L [67] parton distribution functions, the cross sections are calculated in Table 2. 
In this calculation, we neglect the small deviation in sin(/d —a) from unity. Because of the 
small m^, the cross sections of pp —?■ H^A and pp —)■ HA are relatively large as compared 
to those of pp —)■ H^H~ and pp —)■ H^H. We note that the cross section for H^H/H^A 
is about twice larger than that for H~H/H~A, because the parton luminosity of ud in 
the initial proton is larger than that of ud. 

Combining the discussions of the decay and the production of the extra Higgs bosons. 
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we can consider the following processes: 

pp^ H^A/H^H pp^ H^A/H^H + 

pp —)■ HA —)■ 4r, pp —)■ HA —)■ 4r + Z. ^ 

The H^H~ production may not be so useful for the feasibility study of the extra Higgs 
bosons as compared to the above processes because of the small cross section as seen in 
Table 2. The cross sections of 4r 3r (as,-), 4r plus W {a^rw) and 4r plus Z ((J 2 rz) 


can be estimated as follows: 

(jir — O'ha X Bt{H -a- tt), (65) 

03 t ^ (oh+a + Oh-a) X Br(Tr=^ -)■ t ^ p ) 

+ {oh+h + oh-h) X Br ( iJ =^ -)■ t^p) x Bt{H -a- tt), (66) 

OatW - {oh+a + Oh-a) X Br ( i ^'=^ AW^) 

+ {oh+h + Oh-h) X Br(if^ —)■ AW^) x Br(hf —> tt), (67) 

OirZ - Oh A X By{H -a- AZ), (68) 


where we used Br(H -A- tt) ~ 100%. The cross sections of the above processes are also 
shown in Table 2 in the case of tan/? = 35. 

The signal and background for the four and three tan hnal states (without a gauge 
boson) were studied in Ref. [68] in the Type-X 2HDM. It was clarified that the main 
background for these processes can be signihcantly reduced by requiring the high multi¬ 
plicity of charged leptons and tan-jets with appropriate kinematical cuts in the hnal state. 
In this paper, we only show the signal cross sections of the above mentioned processes 
as given in Table 2. Although the detailed background simulation must be necessary to 
clarify the feasibility to detect the signal events, such a study is beyond the scope of this 
paper. 

5.2 Phenomenology of the SM-like Higgs boson 

Another important impact on the Higgs phenomenology in our scenario is found in the 
property of the SM-like Higgs boson h. Because the properties of h; e.g., the width, 
the branching fractions, and the couplings will be precisely measured at future collider 
experiments such as the LHC Run-H, the high luminosity LHC, and the ILC [69], it must 
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mH± [GeV] 

O'H+H- 

(^H+H 

(^H-H 

O'h+A 

O'h-A 

O'AH 

(J 47 - 

^3t 

0'4tW 

0'4tZ 

200 

18.6 

22.0 

11.3 

116 

67.0 

101 

29.3 

50.1 

143 

70.7 

250 

8.0 

9.7 

4.7 

53.5 

29.5 

45.1 

7.2 

12.8 

72.5 

37.4 

300 

3.9 

4.8 

2.3 

28.2 

14.9 

23.2 

2.3 

4.3 

39.4 

20.6 

350 

2.1 

2.6 

1.1 

16.2 

8.2 

13.0 

0.9 

1.7 

22.9 

12.0 


Table 2: Cross sections of the electroweak production processes expressed in Eq. (63), 
and those of the multi-tan processes expressed in Eqs. (65)-(68) at ^/s = 14 TeV in the 
unit of fb. We take = 20 GeV, sin(/l — a) = 1 and tan (3 = 35. 


be quite important to study the deviation in the property of h from the SM prediction. 
In particular, studying the pattern of the deviation in the various h couplings can be a 
powerful tool to determine the structure of the Higgs sector^. 

As we discussed in Sec. 2 , the value of sin(/l —a) describes '^SM-like ness" of h, namely, 
all the h couplings to the SM particles become the same as those in the SM prediction 
in the limit of sin(/l — a) —)■ 1. In other words, once sin(/l — a) 7 ^ 1 is given, both the 
hVV and hf f couplings deviate from those of the SM values. In our scenario, the value 
of sin(/l — a) is determined from Eq. (60). Thus, a small but non-zero deviation from the 
SM-like limit is given. 

In order to describe the deviation in the h couplings, we introduce the so-called scaling 
factors defined as Kx = dhxx/gTxx deviation from unity; i.e, Ak-x = 

Eqs. (23) and (24) and the approximate formulae given in Eqs. (28) and (29), we obtain 


Any ~ 


Ak^ ~ 


tan^/3 

2 

tan^ 13 


mt 


1 -l- 

m 


2m\ 


H± 


m 


H± 


K, ~ -1 - 


ml 
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2m]j± 

2m\ 


m\ 


m 


H± 


m 


H± 


m 


(69) 

(70) 

(71) 


H± 


In the upper panels of Fig. 8 , we show the contour plots for Any and AKg, where 
Akx = Kx — 1, on the m//±-tan (3 plane. In the lower panel, we show the mH± dependence 

®In Ref. [70], the pattern of the deviation was investigated in various extended Higgs sectors; e.g., 
models with isospin singlets, doublets and triplets at the tree level. For example, it was shown that the 
four types of Yukawa interactions in the 2HDM can be well discriminated by measuring the correlation 
between the deviation in hdd and M£ couplings [70]. In addition, it was clarified in Ref. [71] that even if 
we take into account the one-loop corrections to the hff couplings, discrimination of the 2HDMs is still 
valid. 
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Figure 8: Contour plots for Any (upper left) and (upper right) on the mH±-tan/3 
plane, where Akx = i^x — The mH± dependence of is shown in the lower panel 
with tan/? = 35. We take = M = mH± and = 20 GeV in all the panels. The 
horizontal dashed line represents the bound from the signal strength using Eq. (73). 

of Ki instead of showing contour plots, because the tan (3 dependence of Kg, can be neglected 
as seen in Eq. (71). For dehniteness, we take tan/? = 35 in the plot for mH±-Kg. We 
find that the deviations in the hVV and hqq couplings are respectively —0(0.1)% and 
— 0(0.01)% which can also be estimated from Eqs. (69), (70). For the Ml coupling, we 
hnd that its magnitude is maximally about 1.6 times larger than the SM prediction, and 
its sign is opposite to the SM one [21]. From the measurement of the signal strength 
of the h —)■ rr channel, he., /Xt-t at the LHC, the magnitude of k^ is constrained. The 
definition of the signal strength is given as 

_ d/, X Br(h ^ XY) 

[ffixBr(ft^Xr)]s„’ ' ' 

where ah and Br(h — XY) are respectively the production cross section of the SM- 
like Higgs boson h and the decay branching fraction of the h —)■ XY mode. In our 
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parameter set, an is almost the same as that in the SM, because of the small A/tg and 
A/ty, so that the signal strength is simply given as the ratio of the branching fraction as 
/i^y ~ Br(h ^ XY)/Bi{h ^ AF)sm. 

The ATLAS and CMS collaborations report the signal strength as Hrr = 1-43 ^q;37 [72] 
and Htt = 0.91 ± 0.28 [73], respectively. By taking a naive combination of them®, we 
obtain 


= 1.08 ± 0.23. (73) 

Thus, the region with \k(\ > 1.27 is excluded at 2(j level, which corresponds to the 
constraint of ^ 270 GeV as seen in Fig. 8 . 

The deviation in the h couplings makes a difference in the branching fraction of h from 
the SM prediction. In Fig. 9, we show the branching fraction of h as a function of mH±- 
Because only the magnitude of the Mi coupling can be larger than the SM prediction, 
only the branching fraction of the h ^ tt mode is enhanced. On the other hand, that 
of all the other modes shown in this hgure are reduced. The size of deviations becomes 
smaller as mH± increases. We note that the branching fraction of the h —)■ fifi decay is 
also enhanced similar to the h ^ tt mode. When mH± = 300 GeV, and all the other 
parameters are taken to be the same as in Fig. 9, Br(/i —)■ fifi) is about 3.4 (2.5) x 10 “"^ 
in the Type-X 2HDM (SM). 

It is important to mention here that there appears the exotic decay mode h —)■ AZ 
in our parameter set as seen in Fig. 9. Although the coupling constant of the hZA 
interaction is proportional to cos(/d — a) which is suppressed by cot (3 as seen in Eq. (29), 
its branching fraction is not so small, especially for the case with small For example, 
we obtain Br(/i —)■ AZ) ~ 7% with = 10 GeV and mH± = 200 GeV. This new decay 
mode of h gives the additional contribution to the four-lepton channel in the SM Higgs 
boson search if A decays into a pair of muon. In the present scenario, Br(A —)■ nfi) ~ 
Br(A —)■ rr) x (m^/m,-)^ ~ 0.0036 is obtained. On the other hand, from Higgs boson 
searches at the LHG, the upper bound on Br(/i —)■ ZA) x Br(A —)■ ii) with i = e or fi 
is given as 5-10 x 10“"^ for 12 < rriA < 34 GeV [43]. This limit is translated into the 

®To obtain the naive combination, we treat the signal strength from ATLAS as 1.43 ± 0.40 by taking 
the average of errors. 
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Figure 9: The branching fraction of h as a function of in the case of = M = mH± 
and tan(3 = 35. The top, middle and bottom panels show the cases with = 10, 20 
and 30 GeV, respectively. The solid (dotted) curves shows the prediction in the Type-X 
2HDM (SM). 
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Figure 10: Ratio of the branching fraction Br(h —)■ 77 )/Br(h —)■ 77 )sm in our scenario 
with tan/3 = 35. The solid, dashed and dotted curves show the cases with mA=10, 20 
and 30 GeV, respectively. The horizontal dashed line shows the bound from given in 
Eq. (80) at 2cr level. 

bound Br(h —)■ AZ) < 14-28% in the Type-X 2HDM. The typical size of Br(h — AZ) is 
below the upper bound as explained in the above. In addition to this channel, e and fi 
are produced from the leptonic decay of r. Thus, the ZA —)■ ilrr —)■ 4£ -f Ipx channel can 
also contribute to the four lepton channel even though the invariant mass distribution of 
the four lepton system is different from that hy ZZ* —)■ M. This will be a subject of a 
future work. 

Next, we discuss the one-loop induced h —)■ 77 decay mode. Because of the 
contribution, the decay rate can be signihcantly modihed even if the h couplings are not 
changed so much from the SM prediction. We note that the deviation in the h^^ coupling 
can be neglected in the decay rate of the h —)■ 77 mode, because its effect appears in the 
tan loop contribution, but the tan Yukawa coupling is too small as compared to the top 
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Yukawa coupling. The decay rate of the diphoton mode is given as 


r(fc ^ 11) = + ^('}A(„(t,) + ^- 4 „(r„a , 

9 

771 

*■“' = 4 ^’ 

where the Ai, A'^^^ and Aq terms correspond to the W boson, fermion and loop 
contributions, respectively. Each of the loop functions are given by 

Ai{t) = -r“^[2r^ + 3r + 3(2r - l)/(r)], (75) 

^V 2 (^) = 2Y/g5r-2[r + (r - l)/(r)], (76) 

7lo(r) =-r-2[r-/(r)], (77) 


where /(r) is dehned as 


/M = 



< 1 ) 
(r > 1 ) 


(78) 


We note that in the SM-like limit, which gives a good approximation for the numerical 
study in our scenario, we obtain the numerical value of the SM W and top loop contri¬ 
butions; he., Ai{Ty^) + A\i^{t^ ~ —6.45. 

We here discuss the impact of the loop contribution to the h —)■ 77 decay. In 
our scenario, since the light CP-odd Higgs boson is required to explain the muon g — 2 
anomaly, the relatively large mass difference between A and is needed to avoid the 
current experimental constraints, which is generated by the Higgs quartic couplings. This 
can be clearly seen by rewriting the A 3 coupling appearing in Eq. (74) as follows 


= 2(m5,± - m7 + XhAAV^, 


( 79 ) 


where the second term is now zero according to our benchmark scenario shown in Eq. (60). 
In addition, the loop function Aq is expanded as Aq{x) = 1/3 -|- 8a;/45 -|- 0{x^) under 
a; 1 or equivalently mH± 3> nih- Thus, even if is relatively heavy, its contribution 
is not decoupled as we can see the asymptotic behavior A 3 n^/( 2 m^±)Ho(rj:^±) —)■ 1/3 
{mH± 00). 
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The ATLAS and CMS collaborations report the signal strength as fx^-y = 1.17±0.27 [74] 
and and fiyy = 1.12 ± 0.24 [73], respectively. By taking naive combination of them, we 
obtain 


fiyy = 1.14: ±0.18. (80) 

In Fig. 10, we show the ratio of the branching fraction of the h —)■ 77 mode as a 
function of mH± with tan/? = 35. We find that the deviation in the branching fraction 
from the SM prediction is obtained in the range of —30% to —15%. The expected accuracy 
for the measurement of the decay rate of the diphoton mode is around —10% at the LHC 
14 TeV 300 fb“^ and 5% at the ILC [75]. Therefore, our scenario is also probed by 
detecting the deviation in the h —)■ 77 decay rate in addition to the extra Higgs boson 
searches discussed in Sec. 5.1 and the measurement of Ke. By looking at the horizontal 
line representing the bound from fiyy, we see that mH± % 220 GeV is excluded in the case 
of = 20 GeV, which weaker than the constraint of k^. 

6 Conclusion 

In this paper, we have explored the possibility to explain the muon g — 2 anomaly in the 
Type-X 2HDM. We have shown in Fig. 6 that the measurement of the leptonic tan decay 
gives an important constraint on the parameter space. As a result, the region which can 
explain the discrepancy in the muon g — 2 aX the la level is excluded by the constraint 
from the tan decay, and that at the 2a level is allowed. We have found that the parameter 
space with 10 ^ rriA ~ 30 GeV, 200 ^ mH,H± ~ 350 GeV and 30 ^ tan < 50 is favored 
by the explanation of the anomaly at the 2a level. 

After hnding the viable parameter region for the muon g — 2, we have discussed 
the implication of the favored parameter region to the collider phenomenology. In our 
scenario, the 4r, 3r, 4t ± Z and 4r + IF signatures are expected from the electroweak 
productions of the extra Higgs bosons at the LHG. The cross sections of these signals are 
shown in Table 2. 

Finally, we have investigated the possible deviation in the property of the SM-like 
Higgs boson h. We have found that the value of the Mi coupling is predicted to be the 
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SM prediction times about —1.6 to —1.0, and the current data of the signal strength 
at the LHC excludes mH± ^ 270 GeV in the case of mH± = ran and tua = 20 GeV. We 
also have evaluated the branching fraction of the h —)■ 77 mode which is modihed by the 
non-decoupling charged Higgs boson loop effect. We have shown that the deviation in the 
branching fraction of the h —)■ 77 mode is about —30% to —15% which can be detected at 
the future collider experiments such as the high-luminosity LHC and the ILC. Therefore, 
the precise measurements of the branching fractions of h —)■ 77 and h —)■ rr gives an 
important indirect test to probe our scenario. Furthermore, we have found that the 
branching fraction oi h ^ ZA becomes a few percent, and it also would provide another 
important signature to test our scenario. Therefore, in the Type-X 2HDM motivated by 
the explanation of the muon g — 2 anomaly, there are various characteristic deviations in 
the 125 GeV Higgs boson property, which will be tested at collider experiments in the 
near future. 
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A Decay rates for the extra Higgs bosons 

The decay rates of the extra Higgs bosons into the fermion pair are calculated at the tree 
level as 

r(//± ^ //') = V2Gp”^n!\v,a'^x^i\x,„^,x,.„^) 

nn ^ //) = Nix'll , (82) 

r(yl ^ //) = (ij.4, x,Ai . (83) 
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where Xab = iV/ is the color factor, and the two body phase space function \{x, y) 

is given by 


\{x, y) = 1 + x‘^ + y'^ — 2x — 2y — 2xy. 


(84) 


We note that the above formulae are applied to all the types of Yukawa interactions. The 
decay rates into the gauge boson associated modes are calculated at the tree level as 

\/2Gp 3 


V{H ^ AZ) = si 


T{H^ AW^) = 




167r 
y/2Gp 3 

IGtt 


{x^H, , 

{x^rp±,Xj^j^±) . 


(85) 

( 86 ) 


B Electroweak production of a pair of Higgs bosons 


Spin and color-averaged parton level cross sections are given as [77] 

2 / 2 2 
^3/2 (^ 

s-I V s ’ s 


(j(s; ud —)■ A) = 


G^m^ SjS-a. 


d{s]ud —)■ H^H) = 


d{s] qq —)■ HA) = 


d(s; qq —)■ H^H ) = 


727r s 
Gpm^ 
72 tt s 


s — m 


w 


X3/2f 


m]j± mjj 


\ s ' s 


Gl 


F^Z/,2 I 2\^y-a 


187r 


ttq;: 


-{v + a )- 


s — rrir 


2 2 

mir m 


y3/2 / '_ZH_ 'HA \ 


9s 


P VqVH s 
„2 ^2 


■^W^w ^2 


2 r;2.,2 


s s 
s 


(87) 

( 88 ) 

(89) 


^tv'^w \s-m% 


X 


A3/2 




(90) 


where vp = —1/2-1- and Vg and Og are dehned in Eq. (35). Production cross section 
at the hadron collider is calculated as 


a{pp —)■ AB) = 2 




-» AB-, S = TS). 


(91) 


where C is the parton luminosity which is dehned as 

r-9iW92(T/i). 

dr L X 


(92) 
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We note that the parton level cross sections of ud —)■ H~ A and ud —)■ H~ H processes 
are the same as those given in Eqs. (87) and (88), respectively. However, because of 
the difference of parton luminosity of the initial proton, the cross section of H^H/A and 
H~H/A are different in the stage after the convolution of the luminosity function shown 
in Eqs. (91) and (92). 
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